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I.  INTRODUCTION 


Lindquist  and  Simmons'  ’ recently  published  a band  model  formulation 
designed  to  handle  radiance  calculations  for  highly  inhomogeneous  and  non- 
isothermal  optical  paths.  An  important  example  of  a case  for  which  a high 
degree  of  inhomogeneity  pervails  is  the  treatment  of  the  propagation  of 
infraied  radiation  from  aircraft  and  missile  plumes  through  long  atmospheric 
paths.  In  many  cases,  the  primary  radiating  species  of  the  hot  gaseous  plume 
are  F^O  and  CO^,  and  the  primary  atmospheric  absorbers  are  these  same 
species.  The  high  degree  of  line  position  correlation  between  the  hot  gas 
emission  spectrum  and  the  cool  gas  absorption  spectrum  precludes  the 
decoupling  of  the  optical  path  into  a source  region  and  an  absorbing  region 
and  performing  a simple  multiplication  of  the  source  band  radiance  by  the 
atmospheric  band  transmittance  to  arrive  at  the  transmitted  radiance  at  a 
sensor  position.  The  entire  optical  path,  beginning  at  the  sensor  position 
and  continuing  through  the  cool  intervening  atmosphere  and  the  hot  gas  source 
region,  must  be  treated  as  a single  radiating -absorbing  media,  i.e.  , calcula- 
tions must  be  performed  on  a single  highly  inhomogeneous  optical  path. 

The  original  formulation  by  Lindquist  and  Simmons^  ^ considered  the 
case  of  a band  of  randomly  arranged  pre  s sure -broadened  Lorentz  lines  of 
equal  intensity.  The  conditions  for  which  a standard  treatment  of  inhomo- 
geneous paths  with  the  Curtis -Godson  approximation  was  likely  to  fail  were 
determined,  and  the  success  achieved  by  the  new  formulation  for  these  con- 
ditions was  demonstrated.  Here,  this  band  model  formulation  is  extended  to 
a randomly  arrang 'd  band  of  Lorentz  lines  distributed  exponentially  in  line 
strength  and  a randomly  arranged  band  of  Doppler  lines  for  both  constant  and 
exponential  line  strength  distributions. 


An  optical  path  is  called  inhomogeneous  whenever  a gradient  in  species 
concentration,  pressure,  or  temperature  exists  along  the  path. 
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II.  RADIANCE  EQUATIONS 


The  mean  radiance  in  a spectral  interval  A^  centered  on  u for  a general 

(2  3) 

optical  path  extending  from  s = 0 to  s = is  ’ 7 


h(u)  = 


T . d t(  v,  s)  , 

L {v,  s)  1 ds 


where  L (v,  s)  is  the  Planck  radiation  function  evaluated  at  the  temperature 
T(s),  and  t(u,  s)  is  the  mean  transmittance  between  the  origin  (s  = 0)  and  the 
path  position  s.  This  expression  is  exact  to  the  extent  that  L ( u,  s)  is 
sensibly  constant  over  &v.  Practical  radiance  calculations  are  generally 
performed  numerically  by  using  the  finite -difference  form  of  the  equation 

Mf)  = L (u,  s.)  |t  (v , s.)  - 7 (v,  s._j)l 

i=l  L J 

and  band  model  approximations  for  the  transmittance  function  ~{v,  s). 

Only  a single  optically  active  specie  along  the  optical  path  is  considered 
in  this  paper.  If  more  than  one  active  specie  is  present,  the  band  transmit- 
tance function  t (v,  s)  can  be  written  as  a product  of  the  band  transmittances 
of  the  individual  species  provided  that  the  positions  of  lines  in  A^  between 
any  twc  species  are  uncorrelated.  Then,  the  general  band  radiance  equation 
can  be  written 


/°  * ™ t 
L(v)  = - I L (v,  a)  f (v,  s)^T 

II  * 


d T.(  J',  8) 
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where 


_ m _ 

t (V,s)  = [J  t (v,  s) 

i=l  1 


and  m is  the  number  of  active  species. 

Within  the  statistical  band  model  formulation  for  an  array  of  randomly 
arranged  lines,  the  band  transmittance  for  a single  specie  can  be  written 


— , , -W(s)/6 

T ( V,  s)  = e 


(2) 


where  W(s)  is  the  mean  equivalent  width  of  all  the  lines  in  Av,  and  6 is  the 
mean  line  spacing  in  At\  In  order  to  provide  a working  model  for  trans- 
mittance calculations,  W (s)  is  generally  computed  by  assuming  an  approxi- 
mate form  for  the  distribution  of  line  strengths  S for  the  lines  within  Av.  If 
P(S,S)dS  is  the  probability  that  a randomly  selected  line  has  a strength  S in 
dS  when  the  mean  line  strength  in  A^  is  3,  then 

W(  s,  S ) = Z"05  P(S,  S)  W ( s , S)  dS 

Jo 


where  W(s,S)  is  the  equivalent  width  of  a single  isolated  line  of  strength  S. 
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III.  HOMOGENEOUS  PATHS 


The  equivalent  width  functions  W for  isolated  lines  of  both  Lorentz  and 
Doppler  shapes  as  well  as  the  mean  equivalent  width  function  W for  these 
same  line  shapes  and  various  assumed  line  strength  distributions  have  been 
studied  extensively  for  the  case  of  homogeneous  optical  paths.  The  case 

of  interest  to  the  present  work  are: 

1.  Isolated  Lorentz  Line 


wl  = 2itV(xl> 


V l = Lorentz  half-width  at  half-height  , 


f(x)  = xe  " Jg(ix) 


i Jj(ix)j 


J = Bessel  functions  of  the  first  kind  of  order  n , 
n ' 


xL  = Su/2ttvl  , 


u = optical  depth  of  the  path 


2.  Band  of  Constant-Intensity  Lorentz  Lines 


P(S,S)  = 6(S  - S) 


6 “ PLtvXL) 


PL  = 2itYL  ^6e  ’ 


This  6 is  the  dirac  delta  function,  not  thj  mean  line  spacing  parameter. 
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& - effective  mean  line  spacing  in  A' 


- average  Lorentz  line  width  in  A v 


xL  = ku/PL  , 


k = S/6 

e 


3.  Band  of  Exponentially  Distributed  Lorentz  Line! 


P(S,S)  = S '*  e'S/S 


6 ' ~ 0Lfe(xL) 


fe  (x)  = x/^1  + 2x 


4.  Isolated  Doppler  Line 


WD  JlnZ  VDg*xD* 


Y D = Doppler  half-width  at  half-height 


g(x)  = 


exp  -xe 


dz 


n n+1 


V (-1)  x 

Z^(n+i  j'v/rnr 


ln2  Su 

77  "vT 


-■ 


5. 


Band  of  Constant-Intensity  Doppler  Lines 


~r=  0Ds<*D> 


« m 1° 

eD  si  ln2  6* 


= average  Doppler  line  width  in  Av 


XD  = ku/PD 


(7a) 


(7b) 

(7c) 

(7d) 


6.  Band  of  Exponentially  Distributed  Doppler  Lines 


8e(*D> 


_ ^ 

e g(xz)  dz 


E 


i i \n  *»+! 

( - 1)  x 


(x  S I only) 


(8a) 


(8b) 


(8c) 


The  function  f(x)  is  the  Ladenburg-Reiche  function.  Simple  approxi- 
mations to  this  function  have  been  discussed  by  Goldman.  ^ An  accurate 


series  expansion  of  about  seven  terms  ; Jr  large  and  small  x for  both  f(x) 

(7)* 

and  g(x)  has  been  given  by  Rodgers  and  Williams.  The  function  g (x) 

(5)  e 

has  been  studied  by  Malkmus,  who  has  given  a series  solution  for  small  x 


and  an  asymptotic  expansion  for  large  x. 


The  summation  term  in  the  large  x expression  for  f(x)  in  Ref.  (7)  should 
be  b a n. 


I 


The  parameter  6 entering  in  these  expressions  is  the  effective  mean 
line  spacing  parameter  and  is  not  necessarily  the  same  as  the  mean  line 
spacing  parameter  6.  The  form  of  the  parameter  depends  on  the  intensity 
distribution  function  assumed  and  is  introduced  to  yield  the  correct  limit  for 
WL/S  in  the  limit  of  strong  line  absorption. 


IV.  INHOMOGENEOUS  PATHS 


A standard  procedure  for  extending  the  band  model  results  for 
homogeneous  paths  to  inhomogeneous  paths  for  the  case  of  the  Lorentz  line 
shape  is  the  Curtis -Godson  approximation. ( The  net  effect  of  this  pro- 
cedure for  isolated  lines  is  to  replace  the  line  strength  and  line  width 
parameters  S and  V ^ that  appear  in  the  expression  for  WL  with  appropriate 
path-averaged  values  of  these  parameters.  For  a band  of  lines,  path- 
averaged  values  of  k and  0^  are  used  in  the  expression  for  W The  use 
of  path -averaged  values  in  the  expression  for  the  equivalent  width  for  any 
line  shape  is  referred  to  here  as  the  Curtis -Godson  approximation.  The 
path-average  definitions  for  an  optical  path  that  contains  a single  optically 
active  specie  for  an  isolated  line  are 


Se(8)  =U(7T  y*8c(s')  p(s')S(s')ds' 


(9a) 


VeL(s)  ~ u(s)S 


J c(s'  )p(s'  )S(s'  ) yL(s'  )ds\ 


(9b) 


W8)  = u(s)  S ( s ) J c(s')p(s')S(s')yD(s')ds'  , (9c) 


S (s)u(s) 

XT  (s)  = 3 7— r- 

L 2^eL(S) 


(9d) 


r / v = se(8)u(8) 

D ~ V tt  Y„n(s) 


(9e) 
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For  bands  of  lines,  they  are 


l 


ke(s)  = u[s)  f c(s' )P<S' )k(s' ) ds' 

0 

PeL(8)  =^)TTi)  / c(s'>P(s'>k(s')pL(s')ds' 
6 *0 

PeD(8)=^yil7)  / =(»')PC')i?(.')6D(.')d.' 

e 0 


(iOa) 


(10b) 


(10c) 


xL(S)  =ke(S)u(s)/3eL(s)  , (lOd) 

xD(s)  =ke(s)u(s)/3eD(s)  • (10e) 

The  function  u(s)  is  the  optical  depth 

u(s)  = /c<.'  )p(s'  )ds' 

0 

where  p(s)  is  the  total  pressure,  and  c(s)  is  the  mole  fraction  of  the  active 
gas  at  the  path  position  s. 

The  essence  of  the  Curtis -Godson  approximation  is  to  use  path- 
averaged  values  in  the  expressions  for  the  equivalent  widths  W and  W. 
Equivalently,  through  equation  (2),  this  substitution  of  averaged  values  is 
made  directly  in  the  expression  for  the  band  transmittance  function  T{\>,  s). 

As  pointed  out  by  Lindquist  and  Simmons,  ^ however,  it  is  not  r ( v,  s)  that 
appears  directly  in  the  band  radiance  equation  (1),  but  rather  the  derivative 
of  the  transmittance  d T(u,  s)/ds.  They  concluded  that,  if  approximations  have 


I 


to  be  made,  better  accuracy  could  be  obtained  by  making  approximations  to 
the  transmittance  derivative  function  rather  than  to  the  transmittance  function 
itself. 


From  equation  (2),  the  transmittance  derivative  is 

dr(  v,  s)  _ ..  1 dW(  s ) 

dS  ‘ ' T(  ’ ' 6 ~~dT~ 

The  approximation  of  the  transmittance  derivative  thus  depends  on  the  approxi- 
mation of  the  mean  equivalent  width  derivative  dW(s)/ds  and  a numerical 
integration  of  dW(s)/ds  to  obtain  W(s).  Then  T(v,  s)  can  be  obtained  from 
equation  (2)  and,  finally,  d7(  l>,  s)/ds  evaluated  for  use  in  equation  (1).  This 
expression  for  the  transmittance  derivative  is  obtained  simply  by  differentiating 
the  statistical  band  model  expression  for  the  mean  transmittance.  By  starting 
from  first  principles  for  a random  array  of  spectral  lines  (see  Ref.  3 for 
example),  however,  it  can  be  shown  that  the  mean  value  of  the  transmittance 
derivative  in  Av  is  given  by  the  same  expression,  i.e., 

dT(  v,  s)  _ dT(  v,  s) 
ds  ds 

As  for  W(s),  the  approximation  to  dW(s)/ds  can  begin  with  an  examina- 
tion of  the  equivalent  width  derivative  for  isolated  lines.  The  definition  of  the 

(3) 

equivalent  width  of  an  isolated  spectral  ine  for  a general  optical  path  is 


W(s) 


_ 

rB 

/ 1 

-00 

1 - exp 

-1  c(  s'  )p(  s'  )k(  v,  s'  )ds' 

- o y 

(H) 


where  k(  v , s)  is  the  spectral  absorption  coefficient  of  the  line  evaluated  at 
path  position  s.  The  equivalent  width  derivative  is  thus 


® r i 

d^L8.).  = c(s)p(s)  f k(  v,  s)  exp  - f 


c(  s'  )p(  s'  )k(  v,  s'  )ds'  ds 


The  spectral  absorption  coefficient  for  general  line  shapes  depends  on  the  line 
strength  S(s)  and  width  y(s).  If  S(a)  and  \(s)  were  replaced  throughout 
equation  (1?.)  by  the  p c :h-averaged  valuep  Sg(s)  and  y (s),  respectively,  the 
result  would  be  the  Curtis -Godson  approximation.  The  Lindquist -Simmons 
approximation  is  carried  out  by  substituting  the  path -ave raged  values  S (s) 

g 

and  ve(s)  into  only  part  of  the  expression  given  by  equation  (12).  Specifically, 
Sg(  s)  and  y (s)  are  substituted  only  into  the  k(  v,  s'  ) term  that  appears  in  the 
exponential  factor. 


A.  ISOLATED  LINES 


The  case  for  the  Lorentz  line  shape  was  considered  by  Lindquist  and 
Simmons.  The  spectral  absorption  coefficient  for  this  line  shape  is 

S(s)  VL(s)  (13) 


kL(v,  s) 


TT  ] .2  2 

(v  - Vq)  + \L(S) 


where  is  the  line  center.  Substitution  of  equation  (13)  into  (12)  gives 


dW  (s) 
L 

ds 


= c(s)p<s)  f j-5 hZh 

4 I Vs1  + I"-  "01  , 


X exp 


i/ 


c(s'  )p(s'  )S(  s'  )\L(s'  ) 
\l<8'  ) + ( v - uQ)Z 


ds‘ 


dll'-  vQ) 


Substitution  of  yeL<s)  ior  V L(  8 ' ) in  the  exponential  term  and  use  of  the 
definitions  for  Sg  [equation  (9a)]  and  xL  [equation  (9d)]  yields 


The  same  result  occurs  also  if  VeL<s)  is  subsituted  for  \l(b')  in  only  the 
denominator  of  the  exponential  term. 
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( VeL(s)  + U - ^ ) 

By  noting  that  the  integral  is  even  with  respect  to  the  integration  over  v - 

and  making  the  substitution  tan  (0/2)  = (v  - i'J/y  T , the  final  result  can  be 

U 6L 

written  as 


dW  (s) 

— = c(s)p(s)S(s)yL(xL,  pL) 


(14a) 


where 


yLU.  p) 


e -x(  1 + cos  0) 

2 2 

(p  +l)  + (p  - 1)  cos  0 


d0 


(14b) 


and  where  pL  = s)/y  s) , The  parameter  pL  is  the  ratio  of  the  local 
value  of  the  line  width  to  the  path-averaged  value  of  the  line  width  up  to  that 
local  point  and  is  the  urimary  index  of  the  degree  of  inhomogeneity  prevailing 
at  postion  s. 

For  an  isolated  Dopper  line,  the  spectral  absorption  coefficient  is 


kn(  v,  s) 


-M 


S(  8 ) 


-ln2(  v-p-0)2/y  p(s) 


(lb) 


by  following  a procedure  similar  to  that  for  the  Lorentz  line,  the  equivaxenc 
width  derivative  can  be  derived  as 


where 

i 


* 


and 


dW  (s) 

- ds  = c(s)p(s)S(s)yD(xD,  pD) 


(16a) 


(16b) 


j 


Pl=  vD(s)/VeD,5)  • 

B.  BANDS  OF  LINES 

The  derivative  of  the  mean  equivalent  width  is  obtained  from  the 
derivatives  for  isolated  lines  by  applying  the  distributica  functions  P(S,  S) 
given  by  equations  (4a)  and  (5a).  The  application  of  these  distribution  functions 
is  not,  however,  as  clearly  defined  as  for  their  application  to  W for  homo- 
geneous paths.  The  derivative  expressions  of  equations  (14)  and  (16)  depend 
not  only  on  the  local  line  strength  S(s),  but  also  on  the  path -averaged  line 
strength  S^s),  Thus,  both  S(s)  and  S^s)  could  presumably  have  distributions 
about  mean  values  S(s)  and  S^s),  respectively.  Ideally,  the  averaging  should 
be  done  with  a distribution  function  of  the  form  P(S,  S,  Se>  5 ).  One  possioility 
would  be  to  write 


J 


nmifiW 


so  that  independent  distributions  could  be  applied  to  S and  S However 

e ’ 

this  form  implies  a complete  independence  of  the  variables  S and  5 . This 

e 

is  clearly  unrealistic  because  Sg  is  defined  as  the  average  value  of  S.  Father 
than  assuming  the  complete  independence  of  S and  Sg,  th ; present  approach 
is  to  as  sun'.:  the  opposite  extreme,  i.e.  , that  S and  Sg  are  completely  depen- 
dent. Although  this  situation  is  probably  not  fulfilled  either,  it  is  certainly 
much  more  realistic  than  the  assumption  of  independence.  The  assumption 
of  complete  dependence  is  equivalent  to  the  assumption  that,  for  each  line 
in  Av, 


S*(sj  = C(s)Si(s),  (17) 

where  C(s)  is  a constant  independent  of  the  line  index  i.  Since  the  only 
variation  of  S1(s)  along  an  optical  path  is  with  temperature,  this  assumption 
implies  that  all  lines  in  A^  have  the  same  temperature  dependence.  Since 
C(s)  is  independent  of  i,  equation  (17)  also  applies  to  the  average  line  S(s). 
Thus,  S (s)  = C(s)S(s).  The  mean  equivalent  width  derivative  is  then  evaluated 
according  to 


dW(  s,  S) 
ds 


P(S,S) 


dW(  s,  S) 
ds 


dS 


(18) 


where  the  S^  term  in  dW/ds  is  replaced  by  C(s)S.  Note  that  this  quantity 
is  the  mean  value  of  the  equivalent  width  derivative  and  not  the  derivative  of 
the  mean  equivalent  width,  although  it  is  to  be  considered  as  such  when  com- 
puting dT/ds  or  W by  numerical  integration. 

The  result  of  substituting  equation  (14a)  into  (18)  and  performing  the 
integration  over  S for  the  constant -line -intensity  distribution  function 
[equation  ( 4a)]  is 


dWT(s) 


1 


rl  a 


= c( s)p(s)k(s)y  ( x_  , p ) 


(19) 


where  y.  (x,  p)  is  the  same  function  as  given  by  equation  (14b),  xT  = k (s)u(s)/ 
PeL^3^  an<^  Pl  = For  a kanc^  °f  constant -intensity  Doppler  lines, 

a similar  procedure  gives 


1 dWD(s) 


« "~dY~  = c(s)p(s)k(s)yD(xD,  pD) 


where  ’^(x,  p)  is  the  same  function  as  given  by  equation  (16b),  x^  = k (s)u(s)/ 
PeD^s^  and  Pj)  = ^D^S^^eD^S^  For  an  exPonential  distribution  of  line  strengths 
for  Lorentz  lines,  a similar  procedure  in  which  the  intensity  distribution 
function  of  equation  (5a)  is  used  gives 


j <*WL(s) 


- — - = c(s)p(S)k(s)yL(xL,pL) 


(21a) 


yL(*.  p) 


= l£ 


¥f 

**  r\ 


'q  [( + 1 ) + (p^  - 1 ) cos  9]  f 1 + x(  1 + cos  0)]‘ 


(21b) 


For  an  exponential  distribution  of  line  strengths  for  Doppler  lines. 


1 cTWd(s) 


6 ~dT~  = c(s)p(s)k(8)yD(xD,  pD) 


(22a) 


’rD(x’ 1,1  = j;  f t — 

n 1 + 


(22b) 


DISCUSSION 


In  the  Curtis -Godson  approximation,  the  substitution  of  path-averaged 
values  in  the  expression  for  the  equivalent  width,  or  throughout  the  expression 
for  the  derivative  of  the  equivalent  width,  results  in  a complete  separation  of 
the  spectral  and  spatial  integration  operations.  In  all  cases  for  the  Lindquist - 
Simmons  approximation,  these  coordinates  are  not  separated.  This  result  is 
demonstrated  by  the  fact  that,  for  each  new  path  coordinate  s,  y(x,  p)  must 


be  recomputed  by  a 8 or  z integration.  Except  for  the  case  of  y for  an 

L 


exponential  line  strength  distribution,  closed -form  solutions  for  the  y(x,  p) 
functions  have  not  been  found.  Ir  a practical  application  of  the  formulation, 
the  numerical  reevaluation  ot  y(x,  p)  for  each  nsw  path  position  could  con- 
ceivably be  excessively  time  consuming  even  for  computer  application.  The 
approach  taken  in  the  calculation  of  y{x,  p)  for  practical  application  has  been 
to  perform  interpolations  with  respect  to  x and  p on  a table  of  precomputed 
values  of  y(x,  p).  The  equivalent  width  derivative  functions  y jjx,  p),  y^x,  p), 
and  yD(x,  p)  are  presented  in  Tables  1,  2,  and  3,  respectively.  The  functions 
are  also  presented  graphically  in  Figs,  lb,  3b,  and  4b,  respectively,  as 
functions  of  x for  various  values  of  p.  Details  of  the  numerical  integration  of 
these  functions,  simple  approximations  to  these  functions  for  limiting  values 
of  x and  p,  and  the  closed -form  solution  of  y (x,  p)  for  an  exponential  intensity 
distribution  are  presented  in  Section  V. 

Inspection  of  the  general  equation 


= c(s)p(  s)S(  s)y(x,  p) 


forces  the  following  physical  restrictions  to  be  placed  on  the  function  y with 


respect  to  the  parameter  x.  ^ For  small  x,  the  equivalent  width  must  grow 


linearly  with  x,  i.  e.  , as  x — 0,  y — 1.  Since  dW/ds  can  never  increase 
faster  than  it  does  in  the  linear  region,  y must  be  a nonincreasing  function 
with  x (and  hence  y < 1 for  all  x).  Clearly,  y must  also  remain  positive  for 
all  x.  All  of  these  features  are  evident  for  the  y(x,  p)  computed  in  the 
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Table  I.  Equivalent  Width  Derivative  Function  for  an  Isolated  Lorentz  Line 
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Table  2.  Equivalent  Width  Derivative  Function  for  an  Isolated  Doppler  Line 
or  Band  of  Equal- Intensity  Doppler  Lines  (Continued) 


Table  3.  Equivalent  Width  Derivative  Function  for  an  Exponential  Intensity 
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Lindquist -Simmons  approximation  (Figs,  lb,  2b,  3b,  and  4b).  IntheCurtis- 
Godson  approximation,  the  function  yfx,  p)  can  be  derived  as 

y(*.P>  = (2  - p,^i2!  + (p  . u!!U0  _ (23) 

where  h(x)  is  any  one  of  the  equivalent  width  functions  f(x),  f (x),  g(x),  or  g (x) 
depending  on  line  and  distribution  type  considered.  The  functions  y(x,  p)  in  The 
Curtis -Godson  approximation  are  plotted  in  Figs,  la,  2a,  3a,  and  4a  for  com- 
parison with  the  appropriate  functions  computed  in  the  Lindquist-Simmons 
approximation.  It  is  immediately  evident  that,  for  certain  conditions  of  the 
parameters  x and  p,  the  physical  restraints  imposed  on  y(x,  p)  are  exceeded 
in  the  Curtis -Godson  approximation.  The  dominant  failure  for  all  cases  is 
that  y( x,  p)  exceeds  unity  for  various  ranges  of  x whenever  the  inhomogeneity 
factor  is  much  greater  than  ~2.  For  any  value  of  p,  however,  y(x,  p) 
eventually  becomes  less  than  unity  for  large  enough  x.  Thus,  the  primary 
condition  for  failure  of  the  Curtis -Godson  approximation  is  the  condition  of  a 
large  inhomogeneity,  resulting  in  a large  value  for  p,  existing  near  the  be- 
ginning (small  or  medium  x)  of  the  optical  path.  This  is  the  condition  that 
would  prevail,  for  example,  in  calculations  for  l^O  or  CC>2  flame  radiation 
propagation  through  a cool  intervening  atmosphere  in  a band  wing  region. 

If  the  condition  for  the  failure  of  the  Curtis -Godson  approximation  occurs 
at  any  time,  the  condition  may  be  self-perpetuating  because,  although  the 
dimensionless  optical  depth  parameter  x is  a measure  of  depth  into  the  optical 
path  and  generally  increases  with  the  geometric  path  parameter  s,  there  are 
conditions  for  which  x can  decrease  with  increasing  s.  Differentiation  of 
equation  (9d),  for  example,  yields 


dxL(s) 

ds 


c(s)p( s)k(s) 


- (s) 


r 


R 


[2  - PT  (s)] 


X 


(b) 

Figure  1.  Equivalent  Width  Derivative  Function  for  an  Isolated 
Lorentz  Line  or  Band  of  Equal-Intensity  Lorentz 
Lines,  (a)  Curtis- Godson  approximation;  (b) 
Lindquist- Simmons  approximation. 


Equivalent  Width  Derivative  Function  for  an  Expo 
nential  Intensity  Distribution  of  Lorentz  Lines. 

(a)  Curtis- Godson  approximation;  (b)  Lin  .quist- 
Simmons  approximation. 


Figure  3.  Equivalent  Width  Derivative  Function  for  an  Isolated 
Doppler  Line  or  Band  of  Equal- Intensity  Doppler 
Lines,  (a)  Cut tis-Godson  approximation;  (b) 
Lindquist-Simmons  approximation. 


Figure  4.  Equivalent  Width  Derivative  Function  for  an  Expo 
nential  Intensity  Distribution  of  Doppler  Lines. 

(a)  Curtis-Godson  approximation;  (b)  Lindquist- 
Simmons  approximation. 


I Inis,  whenever  p is  greater  than  2,  Xj  decreases  with  increasing  s and  will 
not  increase  to  a value  where  y is  again  less  than  unity.  Another  failure  of 
the  Curti s -Godson  approximation  occurs  for  Doppler  lines  where  y(x,  pi  be  - 
comes  negative  for  certain  values  of  p less  than  unity. 

Certain  other  features  of  the  results  for  the  y(x,  p)  functions  should  be 
discussed.  In  the  Dindquist -Simmons  approximaii'",n,  *he  changeover  from 
the  small  x to  the  large  x form  of  the  function  generally  occurs  much  more 
rapidly  than  in  the  Curti s -Godson  approximation.  This  is  most  evident  for 
yj.Jx.pl  for  small  p (Fig.  1).  For  a given  line  shape  within  the  Lindquist  - 
Simmons  approximation,  the  changeover  from  the  small  x to  the  large  x form 
of  the  function  is  very  much  more  rapid  for  the  constant  line- intensity  distri- 
bution than  for  exponential  distribution.  This  is  in  keeping  with  the  same  phe- 
nomenon observed  for  the  changeover  of  *W  between  the  weak-line  and  strong- 
line regions  and  is  explained  b/  the  allowance  of  weak  lines  by  the  exponential 
distribution  function . ^ ^ ^ 

The  range  of  values  of  p for  which  Tables  1-3  wore  prepared  reflect  that 
the  primary  application  of  the  formulation  is  to  atmospheric  transmittance 
problems  that  involve  hot  gaseous  emission  sources.  The  parameter  p has 
the  value  unity  at  the  beginning  of  any  optical  path  The  variation  of  p along 
the  path  then  depends  on  the  path  type  and  the  variation  of  -y ( s > or  p(s)  along 
the  path.  In  application  to  atmospheric  slant  paths,  two  conditions  prevail 
for  pressure-broadened  Lorentz  lines.  Fo*-  an  optical  path  that  originates  at 
a high  altitude  and  progresses  deeper  into  the  atmosphere,  p^  increases 
gradually  and  generally  approaches  a final  value  of  ~2.  For  an  optical  path 
that  originates  low  in  the  atmosphere  and  progresses  upward  to  low-pressure 
regions,  p^  may  decrease  to  values  as  srnaT  as  10  (for  an  atmosphere  cut- 
off altitude  of  100  km,  for  example).  Nearly  all  of  this  p variation  is  caused 
by  the  pressure  variation  of  the  line  half-width  given  by 

% . v o . / 273 

^L(S)  = vLp(s)Wt (T>  • 
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where  \T  is  the  line -broadening  parameter  for  standard  temperature  and 
pressure.  Very  little  variation  is  due  to  the  temperature  variation  of  \ ( s) 

or  of  6 (s)  because  of  the  relative  i sother mality  of  the  atmosphere.  The 
Doppler  half -width  is 


i 


" 


t , 


\ 
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and  is  pressure  independent.  Thus,  for  a Doppler  line,  p^  will  never  vary 
much  from  unity  for  an  atmospheric  path. 

For  an  optical  pi.th  that  progresses  from  a cool  region  into  a hot  region, 
ft  may  increase  substantially  near  the  boundary  of  the  region  and  result  in 
local  p values  in  excess  of  100.  Here,  the  primary  cause  of  the  increase  in 
p is  the  strong  temperature  dependence  of  6^(  s).  For  a hot  gas,  many  more 
lines  will  generally  be  important  in  a spectral  interval  A v than  for  a cool  gas. 
Consequently,  the  effective  mean  line  spacing  will  be  much  smaller  for  the  hot 
gas,  and  the  ft  value  will  be  much  larger.  The  degree  of  inhomogeneity  en- 
countered in  passing  from  a cool  region  into  a hoi  region  will  be  about  the 
same  for  both  Lorentz  and  Doppler  lines  because  the  major  variation  is  in 

s)  rather  than  in  the  relatively  less  -sensitive  square  -root  variation  of  half- 
width with  temperature.  The  possibility  of  very  small  p values  obtained  by 
passing  from  * hot  region  into  a cool  region  is  formally  allowed,  but  is  not  of 
much  practical  importance  and,  hence,  is  not  considered  here.  This  limitation 
has  the  practical  advantage  of  not  having  to  consider  the  evaluation  of  y(x,  p) 
for  Doppler  lines  for  values  of  p^  much  smaller  than  the  minimum  value 
likely  to  be  obtained  for  strictly  atmospheric  paths,  i.  e.  , about  0.  5. 

The  variation  of  y(x,  p)  with  p also  requires  some  discussion.  The 
solutions  for  all  of  the  y(x,  p)  functions  display  the  feature  that  y(x,  p)  is 
nondecreasing  as  p increases  and  that  y I as  p For  an  isolated  line, 

an  increase  in  p to  a large  value  at  some  point  along  the  path  means  that  the 
local  value  of  the  line  width  is  larger  than  the  width  averaged  over  all  of  the 
previous  path.  If  the  local  width  is  much  wider  than  the  averaged  width,  the 


line  will  continue  to  grow  essentially  as  ,i  new  line  with  only  a small 
perturbation  near  the  line  center  caused  by  the  previous  path.  Thus,  the 
equivalent  width  will  grow  nearly  linearly  with  y near  unity.  For  a band  of 
lines,  the  same  effect  with  line  width  may  occur  or,  as  is  more  likely,  the 
mean  equivalent  width  will  grow  as  a new  band  because  of  the  addition  of  m 
new  weak  lines  at  the  inhomogeneity,  i.e.  , the  increase  in  p is  caused  by  a 
large  decrease  in  6 rather  than  an  increase  in  y . 


V.  SOLUTIONS  FOR  y(x,  p)  FUNCTIONS 


NUMERICAL  INTEGRATION  OF  y.  (x,  o) 


A straightforward  numerical  integration  of  equation  (14b)  was  performed 
by  dividing  the  region  from  0 = 0 to  0 tt  into  N equal  sized  intervals,  inte- 

( 8) 

grating  over  each  interval  w'ith  a 1 2 -point  Wedle  integration  approximation, 

and  allowing  N to  increase  until  the  desired  accuracy  was  achieved  for  each 

5 

fixed  x and  p.  The  data  in  Table  1 are  accurate  to  at  least  1 part  in  10  . 

Vithout  some  form  of  smoothing  transformation,  this  integration  procedure 

would  not  be  practical  for  an  in-program  computer  calculation  of  y(x,  p)  for 

each  new  path  coordinate  s because,  for  certain  extremes  of  x and  p,  the 

integrand  becomes  very  sharply  peaked  about  0 0 or  8 - tr  or  both.  Analysis 

2 2 

of  the  integrand  reveals  that,  for  x >(1  - p ) / 2 p , the  integrand  peaks  around 

2 

0 tt;  for  x < ( 1 - p )/2,  it  peaks  around  0 0;  and  for  x between  these  limits, 

it  peaks  around  both  0 = 0 and  0 = tt.  For  extreme  conditions  of  these  cases, 
the  integrand  becomes  highly  concentrated  in  the  regions  of  the  peaks.  Although 
detrimental  to  straightforward  numerical  integration  techniques,  thisphe- 
nonema  makes  it  possible  to  find  asymptotic  and  simple  approximate  solutions 
for  y(x,  p)  in  these  extremes.  For  example,  the  asymptotic  solution  for  y(x,  p) 
for  large  x is 


y(x,  p)  ~ 


2„  x + 


The  conditions  for  which  this  expression  is  accurate  to  a given  value  are 

summarized  in  Table  4.  These  conditions  were  determined  by  comparison 

with  exact  numerically  determined  values. 

The  limiting  form  for  y(x,  p)  as  p becomes  small  is  y(x,  p)  — exp(  -2x). 

Combining  this  result  with  the  asymptotic  form  yields  the  following  approxi 

. 4 

malion,  which  is  accurate  for  all  x to  at  least  1 percent  for  p i 10  , 
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y(x.  p)  - e'2x  + • (25) 

The  rapid  changeover  from  the  dominance  of  the  first  term  to  the  dominance 
of  the  second  is  evident  in  Fig.  lb.  Note  that,  for  p - 1,  y(x,  p)  is  the  same 
as  for  the  Curtis -Godson  approximation,  or,  equivalently,  y(x,  1)  = df(x)/dx 
by  equation  (23). 

B.  SOLUTION  FOR  yL  (x,  p)  FOR  AN  EXPONENTIAL 
LINE-INTENSITY  DISTRIBUTION 

The  integral  of  equation  (21b)  can  be  evaluated  in  closed  form  by  the 

method  of  residues.  Since  the  integrand  is  even  in  G,  the  limit  of  integration 

i G 

can  be  extended  to  0 = 2ir.  By  introducing  the  complex  variable  z = e , the 
integral  can  be  transformed  (for  p * 1)  into 


y(x,  p) 




itr(p2  - 1) 


7‘f 


2 , 

z dz 


z2  + 2 K + 1 


kP2  * * 


z + 1 


Z2  + 2/i^\z  + l12 


(26) 


where  the  contour  C is  the  unit  circle  about  z = 0.  The  integrand  has  first- 
order  poles  at 


z 


1 


and  second-order  poles  at 


Zj  = - ^[(1  * x)  ” ^ 1 + 2 x ] , 


z = - —[(1  + x)  + \!  1 + 2xj 


(27a, b) 


(27c) 


(27d) 


Only  the  poles  z,  and  z«  lie  within  the  contour  C.  The  residue  at  z = z is 

1 J 1 


1 


<Z|  - *2>u, 


z3),:|zl 


(28a 


The  residue  at  z = z^  can  be  written 


,(3) 


(7. 


2 2 

Zl)  (z3  - z2>  (Z3  - z4} 


3 Z3(Z1  + Z2)(Z3  + Z4>  ' 2(Z3  + z i Z2Z4 ) 


(28b 


Thus,  the  solution  for  y(x,  p)  is 


y(x,  p) 


(29) 


The  solution  is  too  complicated  to  be  written  out  explicitly  in  terms  of 
x and  p. 

When  p - 1 , the  resultant  integral  is 


y(x,  1) 


zdz 


(30) 


Thus,  only  the  second-order  pole  at  z,  is  relevant.  The  resultant  residue  is 


and  the  solution  for  y(x,  p)  in  terms  of  x is 


y(x,  1)  = 


1 + x 


(i  + 2x) 


372 


(32) 


As  a check,  note  that  this  expression  is  the  same  as  for  the  Curtis -Godson 

approximation,  i.e.,  y(x,  1)  = dfe(x)/ds. 

If  the  contour  integration  technique  is  applied  to  the  function  y (x,  p) 

L 

for  an  isolated  Lorentz  line,  it  is  found  that  the  integrand  has  an  essential 

- x / 2 z 

singularity  at  z = 0 from  a term  of  the  form  e . The  residue  for  such  a 

singularity  can  be  obtained  only  from  a Laurent  expansion  of  the  mtegrand 
around  the  singularity.  Unfortunately,  this  expansion  procedure  is  quite 
complicated  and  eventually  yields  the  residue  in  the  form  of  an  infinite  series. 
Lindquist  reports  that  he  has  obtained  a solution  for  this  integral  in  the  form 
of  an  infinite  summation  of  Bessel  functions. 

C.  NUMERICAL  INTEGRATION  OF  yD(x,  p) 


The  integrand  of  the  integral  of  equation  (16b)  has  a single  maximum  at 

2 - 1 / 
z.-O  whenever  x p < 1 and  a single  maximum  at  z = p yli.Vp  ) whenever 

2 U 

xp  >1.  The  integral  was  evaluated  numerically  in  a manner  similar  to  the 

evaluation  of  y (x,  p).  A value  of  z = z was  determined  such  that  the  integral 
L m 

could  be  accurately  written  as 


yU.p)  = ^ 


’/ 


m 


2 

-z  -xe 


2 2 
-P  z 


1 - 


erf(  z )] 
m 


33) 


For  an  absolute  accuracy  less  than  « , z should  satisfy 


m 


z > 
m 


G.  H.  Lindquist,  private  communication. 


-39- 


For  xp  - 1,  the  integral  was  evaluated  by  numerical  integration 
between  z = 0 and  z = z^  with  the  error  function  remainder  added.  For 
x p > 1,  the  integrals  for  the  regions  z = 0 to  z zQ  and  z = zQ  to 
z zm  were  evaluated  separately  and  then  added  along  with  the  remainder 
term.  The  data  in  Table  2 are  accurate  to  1 part  in  10^  or  better.  For 
P 1,  y(x,  p)  = dg(x)/dx. 

A series  solution  for  y(x,  p)  can  be  obtained  by  writing  the  integrand 
as  a product  of  two  exponential  terms,  expanding  the  term  involving  x and 
p into  a power  series,  interchanging  the  order  of  summation  and  integration, 
and  performing  the  resulting  elementary  integration  with  respect  to  z.  The 
result  is 


v(*.e)  -~x*  , . (34) 

n-0  (n  + 1 ) ! \J  1 + np 

This  expression  is  useful  for  calculating  y(x,  p)  for  small  values  of  x. 

The  limiting  form  for  y(x,  p)  as  p becomes  small  is  y(x,  p)  — e'X.  For 
very  large  x and  xp  > 1,  the  first  term  of  the  asymptotic  expansion  for 
y(x,  p)  is 


\fl  1 

y(X’ P)  ~ ? 772  " / T • (35) 

(xp  e)  P yin(xp  ) 

This  expression  is  useful  foi  evaluating  yD(x,  p)  for  x > 100  and  p < 2. 

D.  NUMERICAL  SOLUTION  FOR  7D(x,  p)  FOR  AN 

EXPONENTIAL  LINE-INTENSITY  DISTRIBUTION 

The  consideration  of  the  integral  of  equation  (22b)  is  almost  identical 
to  that  for  yD(x,  p).  The  integrand  has  a single  maximum  at  z 0 for 
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x(2p2  . 1)  < 1 and  a single  maximum  at  z^  = p v ln[x(2p2  - 1)]  for  x(2p2-l>l 
The  integral  Cc  n be  written  as  a definite  integral  between  the  limits  z = 0 and 
z = z plus  an  error  function  remainder  term  to  an  absolute  accuracy  ( for 

■ 

A series  expansion  for  x £ 1 can  be  obtained  as 

y(x.p).£  . (36, 

n-o  V 1 + np 


A useful  approximation  for  p < 0.  5 is 


y(x,  P)  = — ===== 

(1  + x)3/2  J 1 - x(2p2  - 1) 

This  approximation  is  accurate  to  at  least  2 percent  for  all  x.  As  p goes  to 
zero,  y(x,  p) -►  ( 1 + 2x)  . An  asymptotic  solution  (though  not  too  useful)  for 

largo  x is 


1 17  2 <x3/2  1 

y(x,  p)  ~ ^ 2 — ~~ 2 — > ' ==  • <37) 

2p  [x<2p2  - l)]1/p  v In  [ x( 2p2  - 1)] 

Finally,  y(x,  1)  = dge(x)/dx. 

The  tabular  interpolation  for  y(x,  p)  is  facilitated  for  this  case  by 
interpolating  with  respect  to  the  parameter  q - p/(  1 + p),  which  has  limits 
0 and  1,  respectively,  for  p = 0 and  p = ®. 
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